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O ■ Abstract 
p , 

Neutrino masses arising from the spontaneous violation of ungauged lepton- 
number are accompanied by a physical Goldstone boson, generically called Majoron. 
Qj I In the high-density supernova medium the effects of Majoron-emitting neutrino 

decays are important even if they are suppressed in vacuo by small neutrino masses 
and/or small off-diagonal couplings. We reconsider the influence of these decays on 
the neutrino signal of supernovae in the light of recent Super-Kamiokande data on 
solar and atmospheric neutrinos. We find that majoron-neutrino coupling constants 
in the range 3 x 10^^ ^ 5 ^ 2 x 10^^ or g > 3 x 10^^ are excluded by the observation 
of SN 1987A. Then we discuss the potential of Superkamiokande and the Sudbury 
Neutrino Observatory to detect majoron neutrino interactions in the case of a future 
galactic supernova. We find that these experiments could probe majoron neutrino 
interactions with improved sensitivity. 



1 Introduction 



The solar 0] and atmospheric |^ neutrino problems have provided two milestones in 
the search for physics beyond the Standard Model, giving strong evidence for z/e and 
z/^ conversions, respectively. Although flavour changing neutral current interactions 
(FCNC) 0, ^ ^ and/or neutrino decays |^ may play an important role in the inter- 
pretation of the data [0, |^, we concentrate here on oscillations involving very small 
neutrino mass splittings, because they provide the simplest picture |TD|, |TT], . Moreover, 
reconciling the LSND data [|13| and/or hot dark matter suggests the existence of a light 



sterile neutrino Pending the confirmation of LSND results by future experiments, 
we choose however to use in this work just the solar and atmospheric data assuming the 
absence of sterile neutrinos, FCNC and/or of heavier neutrinos. 

Within the simplest extensions of the Standard Model in which neutrino masses are 
introduced by an explicit breaking of lepton-number via the seesaw mechanism there are 
neutral-current-mediated neutrino decays z/' ^ 3z/ 0, but these decays are extremely 
slow for the neutrino masses of interest for us here. If neutrino masses arise from the 
spontaneous violation of ungauged lepton number, the corresponding Goldstone boson 



15| brings in the possibility of new, potentially faster, 2-body neutrino decays g], 

p'^u + J. (1) 

However it is well-known that in the simplest well-motivated majoron models the neu- 
trino mass and coupling matrices are proportional, so that decays are highly suppressed in 
vacuo^ |]16|. As a result, neutrino decays become cosmologically and astrophysically irrel- 



evant in these models. Moreover, for the small masses indicated by solar and atmospheric 
neutrino data, the annihilation channels u' + u' J + J are also negligible. 

The purpose of the present work is to discuss the possible impact of neutrino- majoron 
interactions on the supernova (SN) neutrino signal. Together with the early Universe 
SN are the only site where neutrinos are in thermal equilibrium and so abundant that 
neutrino- neutrino interactions become important. Therefore, SN physics has been one of 
the main tools to derive limits on neutrino parameters in majoron models |]I8|, |T^, ^ ^ 



The arguments used in the literature can be divided mainly into two classes: The 
first uses the fact that majoron interactions violate lepton number (like in u u + J or 
u + u ^ u + u) and, thereby, can reduce the trapped electron lepton number fraction Y^^ 

^This suppression can be avoided in some models due to a judicious choice of the quantum numbers 
H 0. In such case the effects discussed here, ahhough important, are not essentiaL 
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in the SN core ITS, M . Since in several SN simulations it was found that a minimum value 
of ^ 0.38 is needed in order to have a successful bounce shock, the strength of the 
lepton number violating interactions can be correspondingly restricted. Unfortunately, 
the older works |]TB| using this argument did not take into account the derivative nature 
of the coupling of a Goldstone boson calculating the scattering cross sections. Moreover, 
an approximation valid only in vacuum was used for the neutrino decay rates. In ref. |T^ 
the correct medium decay rates were used for the first time, giving the stringent limit 

r/ee < 5 X 10-6. 

The second argument used is a "classical" energy loss argument: the observed neutrino 
signal of SN 1987A should not be shortened too much by additional majoron emission 
[1T9| , pO| , [2T[| . The most comprehensive work along this line is Ref. I^T] which also uses the 



correct derivative coupling of the majoron. The limits found there depend on the value 
of the {B — L) breaking scale. For example the range 5 x 10"'' < (? < 6 x 10"^ is excluded 
for f = 20 GeV while 1 x 10"^ < ^ < 7 x 10"^ is excluded for v = 500 GeV. The tau- 



neutrino masses considered in Ref. are, compared to the values currently discussed, 
rather large, 100 eV < m^^ < 30 MeV. It is therefore of certain interest to extend their 
discussion to lower neutrino masses presently indicated by the solutions to the solar and 
atmospheric neutrino anomalies. 

These astrophysical limits should be confronted with the available laboratory con- 
straints. While there is a stringent limit from jS/S experiments 



Y.9ijUeiUej <3X 10- 



the limits from pion and kaon decays are rather weak, 

E all < 3 X 10-^ and ^ < 2.4 x 10^ 



(2) 



(3) 



Note also that individual couplings Qij could be larger than the limit from [3 (3 experiments 
due to possible cancellations . 

The main purpose of the present work is the study of the impact of neutrino- majoron 
interactions on the observable neutrino signal of a supernova. Since decays of the type 
i^e + J reduce the z/g flux, a limit on the neutrino majoron coupling constants can be 
derived from the observed signal of SN 1987A. In so-doing one must include in the analysis 
the fact that massive neutrinos may oscillate on their way from the SN envelope to the 
detector. We do that and present the excluded regions for the three currently discussed 
solutions of the solar neutrino problem [0, ^ . Then we estimate the discovery potential 
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of new experiments like Superkamiokande (SK) and the Sudbury Neutrino Observatory 
(SNO) in the case of a future galactic supernova. We find that these experiments could 
probe majoron neutrino coupling constants g down to (7 > few x 10~^. 



2 Matter effect on neutrino-majoron interactions 

We consider the simplest class of models in which neutrinos acquire mass from the spon- 
taneous violation of ungauged lepton number ||15[. In this case it is well-known that the 
massless Goldstone boson J - the majoron - couples diagonally to the mass-eigenstate 
neutrinos z/j to a very good approximation |T^. In other words, after rotation from the 
weak basis z/^ through 6q angle(s) the original coupling matrix g^p transforms into 

Qij ~ 5ijgi . (4) 

We denote by uf^""^ the 4-spinor describing the majorana neutrino with mass rrii and 
helicity h = ±1. 

In this section, we briefly review the effect of a thermal background on the neutrino 
majoron interactions. It was first demonstrated by Berezhiani and Vysotsky in Ref. ||2^ 
that the effective mass induced by the interactions of neutrinos with background matter 
can break the proportionality between the mass matrix niij and the coupling matrix 
gij characteristic of the simplest majoron models in vacuo |jl6|. A thermal background 



consists, except in the early Universe, only of particles of the first generation and, therefore 
distinguishes the electron flavour from the other flavours. 



In Ref. p6|, the Lagrangian describing neutrino majoron interaction in a thermal 
background was obtained in a relativistic approximation similar to the usual treatment of 
neutrino oscillations in matter. Later, the authors of Ref. solved this problem without 



using this approximation and confirmed the results of Ref. in the appropriate limiting 
cases. Here we are concerned with SN neutrinos, which have typical energies around 
10 — 25 MeV. Therefore, we will take advantage of the simpler relativistic approximation 
and we will follow closely Ref. p6| . 

The Hamiltonian iJtot = H{x) describing the evolution of neutrinos may be given as 

-f^tot = Hq + H^cd + Hi^t (5) 
= ^0 + E E ^i'^ V.J-'^V]'^^) + g.M'-h^^j''^ J , (6) 



i,j hi,h 
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where the free Hamiltonian Hq describes the propagation in vacuo, i/med describes the 
effects of matter and ifint takes into account the presence of neutrino-majoron interactions 
which may lead to decays. Instead of using the eigenstates of the free Hamiltonian Hq 
as basis for perturbation theory, we will use the eigenstates of Hq + ifmcd- Therefore, as 
a first step the Dirac equation, {Hq + V)i'j'^^ = idty^^^ has to be solved^. In the case of 
ultra-relativistic neutrinos, it is well-known that a vector-like potential does not change 
the helicity of the neutrinos — its only effect is a rotation of the eigenstates oi Hq + V 
with respect to the mass basis |^0|]. Therefore, the Dirac equation simplifies and reduces 
to the standard form known from neutrino oscillations, 

idtuf"^ = {Hlf + U^^Vo^pUl^^^ (7) 

where Hlf (p + mf / {2p))6ij and V^is is the potential matrix in the weak basis 
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The potentials induced by the charged and neutral currents are Vc = \/2hG pUBiYe + Yj,^) 
and Vat = \/2hGFnB [—\Yn + yv^i where Yi = {rii — n-j)/nB and ub is the baryon 
density. Finally, U is the mixing matrix relating mass and weak basis and defined through 
= Uiafa- Diagonalizing H^^^ + UVW gives the medium states 9^'^^ = 
In the case of a two-flavor neutrino system the mixing matrices, U and f/, can be 
parametrized by 6q and 6^^^ respectively. The diagonalization of the Hamiltonian leads 
us to the following expression for the effective mixing angle, 

sin^ 2^0 + (cos 2^0 - ^^^^ 



2^^^^ (p) = . :T Z ...n. > (9) 



where 



^ ''2pGAY, + Y^JnB' ^ ' 



and Aq = ml — ml. The effective masses in the medium are 



' = -hp{2VM + Vc) + + ml) t , (11) 



2V 1 . 2 

where the upper sign is for mi, the lower for m2 and 



A^^) = J(Ao cos 29q + 2hpVc)^ + (Aq sin 29q)^ . (12) 



^The majoron remains massless because the forward scattering amplitude of Goldstone bosons on 
matter vanishes. 
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For typical neutrino energies and densities near the neutrino-spheres, the parameter 



_g Ao lOMeV lO^Og/cm" 



^ = 6.53 X 10"" ^'^ , (13) 

^ 10-3eV2 p (Ye + Y,Jp ^ ^ 

is much smaller than one. This fact will give rise to the following simplification, 0^'^\p) ^ 

—6*0 and 9^ \p) ^ n/2 — 9q, which, as it will be shown later, will allow us to identify 

medium and weak interaction states. 

In the three-flavor neutrino case the mixing matrix U can be parametrized as U = 

U12U1SU23U0, where the matrices Uij = Uij{6ij) perform the rotation in the ij-plane by the 

angle 6ij and Uq includes possible CP- violation effects In the following we will assume 

for simplicity CP conservation and ^13 = 0, the latter motivated both by detailed fits of 



the atmospheric neutrino anomaly, but also by the results of the Chooz experiment |27 



This simplifies the mixing matrix to Vi = Uiai^a = t^i2f^23^'a p8[ and we can make the 
assignment, 612 = 6q and 623 = Oatm- Notice that for light neutrinos near the neutrino 
spheres the condition \Vaa\ ^ f^f/i'^v) ^ 1) holds and, since in the weak basis the 
potential is diagonal, the medium states can be identified with the weak ones up to a 
rotation in the z/^ — subspace. The expressions can be simplified by choosing this 
arbitrary rotation angle to coincide with 623- Then, only one angle, 6'^^\ will be required 
to connect medium and mass eigenstates, and one will be able to recover the two-flavor 
neutrino case. 

Thus the medium states can be identified with weak interaction states according to Table 
|I] and the coupling matrix cjij in the medium basis can be approximated by the one in the 
weak basis, ga/s- 

Above, we have derived the transformation matrix Uij{0^^'^) between mass and medium 
eigenstates. The relation between the corresponding coupling matrices in vacuo Qif and 
in medium cjij is given by 

/ CfCigu + SfSig22 SiCf g22 - CiSf gu ^ 
CiSf g22 - SiCf gn SfSign + CfCig22 (15) 

V ^33 y 

with Cij = cos[6'*^='=)(pjj)] and Sij = sin[6'*^^^(pi j)]. Notice that the coupling constant (733 
will only appear in decays not involving electron neutrinos or anti-neutrinos, so that it 
will not be important for us. 
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medium state 


weak state 


potential 


i 








£23^^/^ + S23Z^T 






-S23t'M + C23'^r 






C23'^M + ■523'^r 












-S2Zyp. + C23Z^r 





Table 1: Medium eigenstates vf^ equivalent weak eigenstates in the limit ^ 1 and their 
potential energy V = m^^^'^ /{2p). 

Within our relativistic approximation, only helicity-flipping neutrino decays — ^ 
vj + J occur. For these decays, the differential decay rate is 



dpf Sn[ pi )\ 2p, 2pf 



(16) 



In the next section, we will apply this formula to supernova. Additional simplifications 
of the limit ^ -C 1 are in this case ml j / {2pij) = Vij and that the coupling constants do 
not depend on pij. Hence, we can integrate Eq. ([TBI) and obtain as total decay rate 

r = ^(^.-^/) (17) 

and as average energy of the final neutrino {Ef) = Ei/3. Obviously, only those decays 
are possible for which Vi — Vf > 0. 

Besides enhancing the rates for neutrino decay, the dispersive effects of the medium 
open also completely new decay channels of the majoron into neutrinos ||3l|. The majoron 



decays J ^ uf + uf have the same matrix elements like the neutrino decays u^^ —^uf + J 



considered above. Their total decay rate is |3T| 



r = ^(-y, - v.) S 



and now those decays are possible for which Vi + Vj < 0. The symmetry factor S is 
5* = 1/2 if the two neutrinos are identical, and S = 1 otherwise. Finally, we note that the 
neutrinos are emitted isotropically. 
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3 Supernova neutrinos and neutrino majoron decays 



Here we first collect the available limits on majoron neutrino coupling constants from the 
observation of the neutrino signal of SN 1987A and comment on their validity. Then, 
we discuss the sensitivity of new experiments like SK or SNO to probe neutrino- majoron 
interactions in the case of a future galactic supernova. 



3.1 Constraints from collapsing phase 

Massive stars become inevitably unstable at the end of their life, when their iron core 
reaches the Chandrasekhar limit. The collapse of the iron core is only intercepted when 
nuclear density, po ~ 3 x 10^^ g/cm^, is reached. At this point, the implosion is turned 
into an explosion: a shock wave forms at the edge of the core and moves outward. The 
strength of this bounce shock and its successful propagation is extremely sensitive to the 
trapped electron lepton fraction Yl^ = Ye + Y^^, attained by the core during its infall. A 



successful SN explosion occurs only if at least 90% of the initial Y^^ is still present p2 
which translates to 11 (^bounce) ^ 0.375 |^ 



We now use this requirement in order to derive a limit on majoron decays. Such 
decays clearly change Y^^ either by two units (z/g — > z/g + J) or by one (z/g + J)- 

Since the allowed change in Yl is small, we can still use the profiles Ye{t) and Yjj^{t) from 
a "standard" SN simulation [Rl]. One can easily check from Table 1 that only the first of 



the above decays takes place, since Y^ + ?>/2Yy^ < 1/2. Hence the deleptonization rate is 
governed by 

rIV 



dt 

We integrate Eq. (|T9|) numerically from to, the time when neutrinos start to become 
trapped (at p(to) ~ 5 x 10^^ g/cm^), to the time of the bounce tbounce- Note that in the 
numerical simulation whose results we are using Yl increases from yL(to) ~ 0.37 up 



to yL(tbouncc) ~ 0.39. Requiring that ^/.(tbouncc) > 0.375 H, ||, i.e. |AFL,dccays| < 0.015 
we obtain 



gee = gii cos^ 00 + 922 sin' < 2 x 10"^ . (20) 



Using the same argument, the limit (^ee ^ 5 x 10 ® was derived in Ref. [p!9[ 



This limit relies on numerical modeling of SN explosions, in particular on the success of 
the explosion in specific models. However, it is probably fair to say that current supernova 
models have generally problems to produce successful explosions. Therefore, Eq. (pO|) 



8 



cannot be viewed as a trustworthy limit as long we do not have a better understanding 
of SN dynamics. 



3.2 Constraints from majoron luminosity 



Numerical computations of the total amount of binding energy Eb released in a su- 
pernova explosion yield, within a plausible range of progenitor star masses and somewhat 
depending on the equation of state used, 

= (1.5 -4.5) X 10^3 g^g_ (^21) 
This range of values is confirmed by likelihood analysis of the observed Ue spectrum of 



SN 1987A under the hypothesis of small mixing of z/g with other neutrino flavours p6 
Therefore, the parameter space of models which give rise to majoron luminosity large 
enough that the observed z/g signal is significantly shortened can be restricted. The most 
comprehensive analysis of SN cooling and majoron emission using this argument was given 
in Ref. |2l|. 

Our following analysis has three main differences compared to Ref. First, we 

are interested in relatively light neutrinos (m^ < eV^ <^ 2pV) as suggested by the sim- 
plest interpretation of data from solar and atmospheric neutrino experiments. Therefore, 
medium effects become important not only for decays like z/ — > z/' + J but also for scat- 
tering processes like J + u J + u. Second, we do not include the process J + J J + J 
in the calculation of the majoron opacity. A discussion of why self-scattering processes 



do not contribute to the opacity was given in Ref. fS^] for the case ofz/-|-z/— >z/-|-z/ 
scattering. Note that trapping due to J+J-^J+J scattering prevented Choi and San- 
tamaria from excluding majoron-neutrino couplings for the case of light neutrino masses, 
nii,^ < 100 eV. Third, we calculate the majoron luminosity in the trapping regime in a 
different way. 

Let us discuss first the mean free path of majorons. The main source of opacity for 
majorons are the processes J +J z^ + z/, J + v J + z/, and J + v v (cf. |21|]). Taking 



into account the effective mass of the neutrinos, the corresponding mean free paths inside 
the SN core with radius ro ~ 10 km and density p ~ pnuc are given by. 



/- V + i - . + W = 1.9 X 10 V \— j j c,„- (22) 

r'U + .^. + .) ^ 1.5xl0V(^)^(^)%.n-' (23) 



9 



l-\u + J^u) = 1 X lOV cm-^ (24) 

^ ' ^ V25MeV/ ^ ' 

For the effective neutrino mass inside the SN core, m^g ~ '2pV, we use meg ~ 20 keV. 

Requiring that the mean free path of the majorons is > 10 km, we obtain that they 

escape freely for (7 < 3 x 10^^. 

Next, we consider the luminosity for the case that majorons are not trapped. Then 

they are emitted by the whole core volume with luminosity 

C{u + u^J + J) = 5.5 X 10 V ( ^] erg/s (25) 

C{u^u + J) = 4.8 X lO^^g^ erg/s . (26) 

We require now that the majoron luminosity during 10 s does not exceed the maximal 
theoretical value of the binding energy, Cj < 5 x lO^^erg/s, and obtain the limit < 3 x 
10~^. In the trapping regime, the energy loss argument therefore excludes the band 

3 X 10-' <g<Sx 10-^ (27) 

for vacuum neutrino masses m <^ rrics ~ 20 keV. 

Let us now discuss the case that majorons are trapped, g > 3 x 10~^. We recall first 
the standard treatment as presented, e.g., in Ref. |21|. The main assumption is that Cj 
can be approximated by blackbody surface emission 

Cj = —R^jT^ (28) 

of a thermal majoron-sphere which radius Rj is defined to be at optical depth 2/3, 

drl-^{r) = -. (29) 
R.J 3 

Assuming furthermore a temperature profile T(r) outside the supernova core, one can 
determine for which range of coupling constants Cj exceeds a certain critical value. 

There are two shortcomings in this argumentation. First, the profile T(r) oc 
used for r > tq in I^T] represents the temperature of ordinary matter (nucleons, e^, 
photons). However, majorons couple mainly to neutrinos and would at best be in thermal 



equilibrium with neutrinos, if at all, but not with nucleons. Second, the use of Eqs. ( p2| - 
p4|) which were derived for isotropic distributions inside the core is not justified for the 
highly non-isotropic distributions outside the neutrino- or majoron-spheres. Using for 
the density ni{r) of a particle species i with sphere radius Ri the expression 
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and introducing the additional factor (1 — cos 6) ~ {R-JrY which represents the averaging 
over the angle 9 between the momenta of the radially outgoing test majoron and the 
particle i in the mean free path one finds that ?2j(r) oc This renders support to 

our claim that there is a sharp drop in the majoron opacity in crossing the neutrino 
spheres. Note also that it was recently stressed in Ref. [^] that a naive application of the 
Stefan-Boltzmann law can be dangerous for the calculation of neutrino-sphere radii. 

We prefer therefore not to rely on the argumentation presented above. Instead, we 
calculate the volume rather than the surface luminosity, but taking into account that only 



majorons emitted from the shell [tq — I : tq] can escape [|T9|. We find that the majoron 
luminosity drops below C^^tot ~ -Ef,/10 s for (7 > 2 x 10~^. Going even to larger values 
of g, majorons and neutrinos are becoming so strongly coupled that it is reasonable to 
assume that the SN core emits roughly the same luminosity in form of majorons as in 
one neutrino species, ie. Cj ~ £,^^tot/6. Therefore, majoron emission in this regime is 
not constrained at all by the energy loss argument. Combining the both limits obtained, 
majoron couplings one concludes that the range 

3 X 10-^ < ^7 < 2 X 10-^ (31) 

is excluded for vacuum neutrino masses m <C ~ 20 keV. 

A final remark is in order. We have been using rather loosely only one coupling 
constant g in this section. More precisely, we mean by g the element of the coupling 
matrix g^p in the weak basis with the largest absolute value. 

3.3 Constraints from neutrino spectra 

During the Kelvin-Helmholtz cooling phase, t ^ 1 — 10 s after core bounce, the protoneu- 
tron star slowly contracts and cools by neutrino emission. The neutrino luminosities are 
governed by the energy loss of the core and are, therefore, approximately equal for each 
type of (anti-) neutrinos. Since the opacity of, e.g., z/^^^ is smaller than of z/g, due to their 
smaller cross section, their energy-exchanging reactions already freeze out in the denser 
part of the protoneutron star. Hence, one expects the spectral temperatures of Ve to be 
smaller than the one of = {i^^,r, ^'/^.t}- Typically, the average energies (Ei) found in 
simulations [0^ are 



{E^J ^ 11 MeV , (EpJ ^ 16 MeV , (^^J ^ 25 MeV . (32) 
It is convenient to define three energy spheres RE,ue: RE,ue and RE,Vhi outside of which 
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only energy- conserving reactions contribute to the neutrino opacity^. These reactions, 
mainly neutral current neutrino-nucleon and neutrino-nucleus scattering, do not change 
the neutrino spectra, although the neutrinos still undergo several scattering as they diffuse 
outward. Eventually, the neutrinos reach the transport sphere at Rt,ui, where also the 
energy-conserving reactions freeze out, and escape freely. This sphere, i.e. the surface of 
last scattering, is what most authors mean with "the" neutrino-sphere. 

For the numerical evaluation of the decay rates, we use profiles for p(r) and Ye{r) from 
Wilson's SN model at the time t = 6 s after core bounce. For all radii of interest, 
only anti-neutrinos can decay into neutrinos. Correspondingly, the allowed majoron decay 
channels can be determined from Table 1. For the average position of the various neutrino- 
spheres, we use the implicit definition 

piRE,uJ = 2 X 10i3g/cm' , piRE,pJ = 2 X 10i2g/cni3 ^ p{Re,uJ = 2 x 10"g/cm=^ , 

(33) 

and 

p{Rt,uJ = p{Rt,,J = p{Rt,uJ = 2 X 10"g/cm^ . (34) 

In calculating the effect of majoron-emitting neutrino decays on the emitted u spectra 
we must distinguish three regions. If a Ui is produced inside its own energy sphere, r < 
Re,p,, it is still in thermal and chemical equilibrium with the stellar medium. Therefore, 
the production (or the decay) of a Ui within r < Re,u, does not influence the emitted z/j 
spectrum. In contrast, the spectrum of z/j neutrinos is fixed for r > RE,Ui and both its 
decay or its production changes its spectra. The only difference between the two regions 
RE,vi < r < Rt^p- and r > Rt^p. is the different "effective" velocity v of the decaying 
neutrino. In the latter, the neutrino can escape freely (f = 1), while in the first it diffuses 
outward with v ^ X/ {RE,Ui ~ Rt,ui)- Here, A is the average mean free path of the neutrino 
Ui. Hence, we can compute the survival probability N{iJi) of a Ui neutrino emitted from 
its energy sphere as 

{rRt.a, rjr' roo 
- —T,Xr')- rfr'rp,(r') , (35) 

where 

r., = E Tiui-^ui + J) (36) 

•^Neutrino-matter interactions are energy-dependent and, therefore, the concept of a neutrino-sphere 
is evidently an over-simpHfication. 
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is its total decay rate. Denoting with Pij[ri, the probabihty that the decay Ui —^uj + J 
happens in between ri and r2, it follows that 

P^i[ri, ra] = / — r(z7, z/,- + J; r')N{r') . (37) 

and J2j Pij[i^ijf'2] + = 1, where N denotes the neutrino survival probability. 

Finally, we have to consider the effect of majoron decays J ^ ui + uj on the neutrino 
spectra, li g <3 x 10^®, majorons escape freely. Thus they are emitted by the core with 
luminosity given by Eq. (|25|). For their mean energy we assume (Ej) = 3T ~ 50 MeV, so 
that the number of emitted majorons per unit time is Nj = Cj/ (Ej). In the intermediate 
regime, 3 x 10^^ ^ fl' ^ 2 x 10^^, we use the same average energy for the majorons but 
use the reduced luminosity due to emission from the shell [vq — I : vq]. Finally, for 
(7 > 2 X 10~^, we identify simply the majoron-sphere Rj with the energy sphere of the 
neutrino species which has the largest coupling to the majoron. In the following, we will 
use Rj = RE,Uh we will also assume that its average energy is similar to the one of 

{Ej)^{E,J^ 25 MeV. 

The decay probability P = 1 — N oi a majoron is given by 

N{J) = exp |- dr' rj(r')| (38) 

Together with the majoron luminosity Cj this allows us to calculate the spectra of the 
produced neutrinos. 



3.3.1 Neutrino signal from SN1987A 



The existing observational data on the supernova SN 1987A were already used in Ref. |^ 
in order to constrain the allowed permutation between the z/g and all other types of 
neutrinos as a result of neutrino oscillations. At 99% (95%) confidence level, it was found 
that no more than 35% (23%) of the Ue flux could be converted into Uf^^r or u^^r- This limit 
can be applied to any process which diminishes the flux of Ue and is only weakly energy 
dependent, such as the case of medium-induced neutrino decays. Thus one can apply 
their results in order to restrict the neutrino decay models under consideration here. 

Since and u^- feel the same potential, the decay rates r(z/e i^fj.+J) and r(z/e 
J) differ only due to their different coupling constants. We define therefore g'^f^ = g1^ + g1^ 
and present the region excluded by the SN 1987A signal in Fig. |I| in the gee — geh plane. 
The region excluded at 95% confidence level corresponds approximately to the one given 
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by the condition 

^?ee + ^7eM + ^?er>lxlO-^ (39) 

Up to now, we have neglected that massive neutrinos do not simply decay but may 
also oscillate on their way. In the high-density region where the decays occur this ap- 
proximation is legitimate because the medium states coincide essentially with the weak 
flavour eigenstates. 

However, the SN 1987A neutrinos must also propagate first through the SN envelope, 
then through vacuum before they cross the Earth on the way to the detectors. Let us con- 
sider now what will be the impact of neutrino oscillations for the three popular solutions 
the solar neutrino problem, namely small-angle (SMA) MSW, large-angle (LMA) MSW 
and the just-so or vacuum oscillations^. These solutions are characterized by particu- 
lar values of Aq and sin^ 26*0 , which define different regimes in the dynamics of neutrino 
propagation through the supernova. 

Let us first analyze the two MSW solutions. For the usual neutrino mass hierarchy, 
anti-neutrinos will not encounter an MSW resonance on their way through the SN enve- 
lope. Furthermore, they propagate adiabatically, because in both cases Aq ~ 10~^ eV^. 
Therefore, neutrinos that had been created as i^i{p) ~ i^e leave the star as I'lip = 0) = Ui, 
i.e. as a definite mass eigenstate. Consequently, no neutrino oscillations take place on the 
way from the SN to the Earth and the probability of z/g i^fj. transitions (permutation 
factor) equals Pose = K'^/il^i^)!^ = sin^^o, as long as matter effects in the Earth can be 
neglected. Since decays z/^ — > z/^ and oscillations z/q are decoupled, the total survival 
probability can be written as 

N = A^decayiVosc = N^ecyil ' ^osc) • (40) 

where Pqsc is the neutrino conversion probability due to oscillations. 

In the SMA-MSW case, sin^ 29o ^ 7 x 10^^, oscillations can be neglected both in the 
SN envelope and in the Earth, so that Pose ~ 0. Thus, one can use directly the results 
obtained above, shown in Fig. |l|. 

In contrast, the matter effect inside the Earth has to be taken into account in the 
LMA-MSW case, for which sin^ 26'o ~ 0.6 []TU|. The permutation probabilities due to 



oscillations is given by. 

Pose = sin 29q — sin 26 sin(2^o - 26) sm\nd/losc) ■ (41) 



We follow closely the discussion given in Refs. [|6[ 
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The distance d traveled inside the Earth by the neutrinos and the average density p are 
different for Kamiokande and 1MB detectors. For Kamiokande we have d = 3900 km and 
p = 3.4 g/cm^ while for 1MB d = 8400 km and p = 4.6 g/cm^. As an approximation 
we will use therefore the average value Pose = (12Posc,Kam + 8Posc,ibm)/20 according to 
the number of events detected in each detector. In Fig. 0, we show the excluded region 
for sin^ 2^0 = 0.6, Aq = 10"^ eV^, and in addition the experimental limit Eq. (^. The 
coupling g22 is fixed by 

g22 = 9iJl + ^. (42) 
V mf 

Therefore, for mi/Ao 0, g22 becomes larger for constant gu. Hence, the limit on gn 
becomes stronger. Note also that we can not fix g^s by the solar neutrino data. Therefore, 
we have set conservatively (733 = 0. As can be seen from the Fig. ^, for the range of masses 
involved, the values gu > 10~^ are excluded at 95% CL. 

The remaining case to consider is the 'just-so' solution, in which Aq ~ 10^^*^ eV^. For 
such small values of the mass splitting, the neutrino propagation out of the SN is non- 
adiabatic. Therefore, neutrinos leave the SN envelope as flavor eigenstates which then 
oscillate on their way to the Earth. Taking into account that in this case the Earth effect 
is unimportant, their averaged permutation probability due to oscillation is simply 

Po,, = ^ sin^ 29o . (43) 



Let us consider now the typical value sin 2^o ~ 0-9 |TT1]. Then, according to the analysis 



of Ref. [^, this case is already disfavored assuming only oscillation. In any case for 
completeness we have plotted in Fig. ^ the regions corresponding to P = 0.55, 0.7 and the 
experimental limit Eq. (0). 

To summarize, the limits we have obtained in this subsection from the observed z/g 
signal of SN 1987A are more than one order of magnitude stronger than the limit from 
pion decay. They require as theoretical input some knowledge about the spectral shape 
of the emitted neutrino fiuences (e.g. {Ey.)) and are therefore to a certain extent model 
dependent. However this model dependence is much weaker than that implicit in the limit 
from the collapsing phase 



3.3.2 Constraints from a future galactic supernova 

The neutrino signal from SN 1987A observed by Kamiokande and IBM confirmed the 
general astrophysical picture of a supernova explosion. However the small number of 
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observed neutrinos prevents a sensitive probe of neutrino properties or of the equation of 
state of super-dense matter. Meanwhile, several new experiments like Super-Kamiokande 
or the Sudbury Neutrino Observatory (SNO) have been constructed and it seems therefore 
interesting to estimate their sensitivity to new physics, such as the majoron-emitting 
neutrino decays, through their ability to determine the spectra of SN neutrinos. 

All detectors in the near future can detect SN only in our own Galaxy and in the 
Large and Small Magellanic Cloud, with however a much smaller rate than expected for 
the Milky Way. Therefore, we use in the following 10 kpc as the reference distance to the 
hypothetical Supernova. Then, according to Ref . |Q , the number of Ue + P ^ n + e'^ 
events expected in SK is 5310, compared to 220 events in all other channels. While a light- 
water detector like SK acts like a Ue filter, a heavy-water detector has a broad response 
to all types of neutrinos: from the total of 781 events expected at SNO, 200 events are 
neutral-current reactions of {=i^fj,, or v.j) with deuterium. In addition, there are 

82 charged-current reactions of v^. with deuterium. The combined spectral information 
of these reactions should contain the major part of the information extractable from all 
data. In the following, we will consider therefore only the reaction z/g + p ^ "n, + e"*" for 
SK, and Vh^d^Vh^n^-p and z/e + (i — p + p + e~ for SNO. 

The predicted number of events has a rather strong dependence on how the cooling of 



the SN is modeled. In Ref. e.g., the expected number of f/i + rf — > z//i + n+j9 reactions 
is 485 in SNO. An important question is therefore how well different astrophysical and/or 
particle physics processes (e.g. neutrino oscillations versus decays) can be disentangled 
using observational data. We will comment shortly on this question at the end of this 
section. 

The neutrino spectrum n^iEy) emitted by the SN can be be inferred in a rather indirect 
way from the experimentally measured energies Ei of the observed events. Let us consider 
first the reaction z/g + p — > "n. + e"*" and z/e + c?— >p + p + e~. Experimentally, the neutrino 
energy Ey is reconstructed in both reactions from the number of photo-multipliers that 
have been triggered by Cherenkov photons emitted by the e^. The probability that a 
certain energy Ei is ascribed to an event i where the has the energy E^. depends in 
principle on the geometry and the details of the detector. For our purposes it is enough to 
take into account the Poissonian nature of the emission and detection of the Cherenkov 
photons. In this case the probability distribution function that the energy Ei is ascribed 
to an with energy E^. is a Gaussian distribution, 

G{E,^ E,) = -=J^ exp (44) 
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with a{Ee) = \^E^E^, where E^r is the energy resolution. Following Refs. [HI], we use 
E„ = 0.22 MeV for SK and E^ = 0.20 MeV for SNO. Apart from its resolution E^, the 
other main feature of a detector is its efficiency e{E). For both experiments, we use the 
approximation e{E) = eQ6{E—Ecnt) with eo = 1 and a conservative value for E^nt = 7 MeV 
for SK and E^nt = 5 MeV for SNO. Combining the Gaussian kernel G{Ei, E^) with the 
detector efficiency e{E), the relation between the energy Eg, and the measured energy 
Ei is 

n{Ei) = eo / dE G{Eu Ee)ne{Ee) . (45) 

In contrast to the two reactions considered above, the neutral-current reaction Vh + d ^ 
Uh + n+p used to detect Uh in SNO does not allow to measure the energy E^. Nevertheless, 
it is possible to reconstruct partially the energy spectrum comparing the observed signal 
in different detector materials with different energy responses (cf., e.g.. Fig. 3 of Ref. [0). 

As next ingredient, we need the relation between the time-integrated neutrino spectra 
ni{E) and the spectra of the secondary e^'s. For the reaction Ue+P n + e~^, the positron 
spectrum is given by 

N 

MEe) = ^^^ap^+p^n+e+(^ + f^)nuAE + /i) (46) 

where Np is the number of target protons in SK and D the distance to the supernova. 
For the cross section, we use 

a^^+p^n+e+ = ooE - -j M - — + I (47) 

with ctq = 2.3 X 10~^^ cm^, where /i = 1.227 MeV is the proton-neutron mass difference 
and TTif. the electron mass. In the case of the two neutrino-deuterium reactions, we use the 
cross-sections tabulated in Ref. and replace Np by the number of deuterium nuclei. 

Numerical simulations of the neutrino transport show that the instantaneous neutrino 
spectra can be described as Fermi-Dirac distributions with an effective degeneracy pa- 
rameters rji The instantaneous neutrino spectra found are pinched, i.e. their 
low- and high-energy parts are suppressed relative to a Maxwell-Boltzmann distribution. 
Janka and Hillebrandt found that during the cooling process the effect of pinching be- 
comes less important. Moreover, the pinching of the instantaneous neutrino spectra is 
compensated by the superposition of different spectra with decreasing temperatures. 

In the following, we assume for the time-integrated energy spectra of the different 
neutrino types 
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with {Ey^) = 11 MeV, {E^J = 16 MeV and {E^J = 25 MeV. For the degeneracy pa- 
rameter, we use the values of the lower end of the range found by Janka and Hillebrandt, 
namely t] = 3 for z/e, ?7 = 2 for z/g and t] = for Uh. Then the relation between the effective 
temperature Tj is given by 3T ^ 0.751{E,yJ for u^, 3T ^ 0.832(£'pJ for and 3T = (-Ej^^^) 
for Uh. 

We have performed a simulation of the expected neutrino signals in SK and SNO in the 
case of a galactic Supernova. Here have assumed gn ~ (722 ~ 933, as suggested by a scheme 



with all three neutrino masses are quasi-degenerate, in the eV range |48|. Moreover, here 



we have chosen the solar and atmospheric mixing angles to be maximal [^. The resulting 
average neutrino signals disregarding neutrino oscillations are shown for several values of 
gn in Fig. ^ and ^ As one could expect from the discussion in Section 3.3.1, for the chosen 
values of g neutrino decays do not influence drastically the Ue signal (Fig. H). In contrast 
the i>h signal shown in Fig. |^ shows two main features. First, in the presence of majoron 
interactions, the spectra show a surplus of low-energy {E < 18 MeV) and a deficit of 
high-energy {E > 18 MeV) Uh compared to the reference Standard Model spectrum with 
g = 0. One sees that the decays Uh i^i + J reduce the Uh fluence by at most a factor 
2. Since the decays Vh vi + J and the majoron decays J u + u' both produce also 
low-energy h'^-,^J., the signal shows a surplus of events at low energies. 

We have not shown the z/g spectra because they are basically unchanged with respect 
to the reference spectrum with g = 0. Since the Ue energy sphere at RE,ue is located at a 
much smaller density than those of z/g and Uh, most (~ 95%) z/g and Uh decay inside Re,!/^ 
where the u^. are still in chemical equilibrium. Therefore, its spectrum is hardly affected 
by z/g and Uh decays. 

Finally, we want to comment on the dependence of our results on the assumed astro- 
physical parameters. First note that, an overall suppression of all three neutrino signals 
could be explained more naturally by an astrophysical reason than by neutrino decays. 
Similarly if a suppression with respect to the expectations of a given SN model were ob- 
served only in the Uh signal, this would indicate an astrophysical explanation, since the 
Uh temperature has the largest uncertainty. It is therefore of importance that part of the 
spectral information of the Uh signal can be recovered by comparing the signal in different 
detector materials. The signature for majoron neutrino decays would be two-fold. First, 
a reduction of the observed Uh fluence compared to the one expected from the observed 
Uh temperature. Second, a non-thermal Uh spectrum with a surplus of low-energy and a 
deficit of high-energy Uh. 



18 



A more quantitative study of the restrictions on the majoron-neutrino couphngs at- 
tainable at SK and SNO would require a detailed likelihood analysis. However, from Fig. |^ 
one expects that the sensitivity of these experiments will be restricted to g < 5 x 10~^. 
Therefore, we conclude that although these experiments could narrow down the allowed 
window for majoron-neutrino interactions, but not close it completely. 



4 Conclusions 

We have reconsidered the influence of majoron neutrino decays on the neutrino signal 
of supernovae in the full range of allowed neutrino masses, in the light of recent Super- 
Kamiokande data on solar and atmospheric neutrinos. In the high-density supernova 
medium the effects of Majoron-emitting neutrino decays are important even if they are 
suppressed in vacuo by small neutrino masses and/or off-diagonal couplings. In contrast 
to previous works, we have considered scattering and decay processes, taking into account 
medium effects for both kinds of processes. The dispersive effects of the dense SN core 
are particularly important since the currently favoured interpretation of the solar and 
atmospheric neutrino data points towards light neutrinos. 

The observation of SN 1987A excludes two parts of the possible range of neutrino 
majoron coupling constants. In the range 3 x 10^^ ^ ^ 2 x 10^^, where g is the largest 
element of the coupling matrix g^js, the supernova looses to much energy into majorons, 
thereby shortening the neutrino signal too much. For larger couplings, the fiuence of 
escaping z/g is reduced due to decays Ue ^ J^i + J ■ Depending on which is the correct 
solution to the solar neutrino problem, different values of g^^ can be excluded: g1^ + 

> 1 X 10~^ for SMA-MSW (Fig. |ip and c/n > 1 x 10"^ for LMA-MSW (Fig. 0). In the 
case of vacuum oscillations, the predicted number of z/g events is already disfavored (even 
in the absence of decays) by the observed number of events and by the current theoretical 
understanding of supernova explosions. The corresponding sensitivities are displayed in 
Fig-i 

Finally, we have discussed the potential of Superkamiokande and the Sudbury Neu- 
trino Observatory to detect majoron neutrino interactions in the case of a future galactic 
supernova. We have found that although these experiments could narrow down the al- 
lowed window for majoron-neutrino interactions, but not close it completely, reaching a 
sensitivity at the few x 10~^ level, as seen from Fig. ^ 
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Figure 1: SN 1987A constraint on the majoron-neutrino effective coupling constants in 
the Qee — Qeh plauc. Parameters corresponding to the SMA-MSW solution to the solar 
neutrino problem are assumed. 
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Figure 2: SN 1987A constraint plotted in the rrii — gn plane. Parameters corresponding 
to the LMA-MSW solution to the solar neutrino problem are assumed, sin^ 29 = 0.6 and 
Ao = 10-^ eV^. 
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Figure 3: Regions of total transition probability plotted in the mi — gn plane. Here 
we assume the just-so solution to the solar neutrino problem, sin^2^ = 0.9 and Aq = 
10~^° eV^ and include both the effects of neutrino decays and oscillation. 
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Figure 4: Signal expected in Super-Kamiokande due to the reaction Pg + p — > n + e"*" for 
a galactic supernova. Here N^, denotes the number of events for g-n — (sohd line) and 
gii = 10~^ (dashed hne) with gu = c/22 = S'ss- 
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Figure 5: Signal expected at SNO due to the reaction 1/^ + D ^ 1/^+ p + n for & galactic 
supernova. Here N^, denotes the number of events for the indicated values of gu with 

911 — 922 — 933- 



28 



